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We present a measurement of the branching fraction of the decay B^ — > XcoK^ from a sam- 
ple of 89 miUion BB pairs collected by the BAR4R detector at the PEP-II asymmetric-energy 
B Factory at SLAC. The Xco meson is reconstructed through its two-body decays to 7r"'"7r~ and 
K+K-. We measure B{B^ -» Xc»K^)yB{xc» tt+tt") = (1.32i;5;27(stat) ± 0.09(syst)) x 10"® 
and B{B^ XcoK^)xB{xcO ^ K+R-) = (1.49i^ j^(stat) ± O.ll(syst)) x 10'^ Using the known 
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values for the XcO decays branching fractions, we combine these resuhs to obtain 
B{B^ ~* XcoK^) = (2.7 ± 0.7) x 10"*. 

PACS numbers: 13.25.Gv, 13.25.Hw 



In the simplest approximation, weak decays like 
B J/iJj K arise from the quark- level process b ccs 
through a current-current interaction that can be writ- 
ten [c7^(l — 75)c][s7^(l — 75)6]. The colorless current 
C7^(l — 75)0, which can create the J/ijj , can create as well 
the P-wave state Xci- It cannot, however, create XcO, Xc2 
or hci so their appearance would have to be ascribed to 
more complex mechanisms. The b — s- ccs process also oc- 
curs through the interaction of two color-octet currents 
jK-cc)jgb) ^ [c(Aj2)7,,(l-75)c][^(Aj2)7ni-75)6], 
where Xa are color SU{3) matrices. The current Jg 
can create a color-octet cc pair in an S state, which can 
then radiate a soft gluon to produce a P-wave bound 
state 0,0- Alternatively, the XcO, Xc2, he states might 
arise from final state interactions that mix the {cc)K 
channel with channels like D^*^Di*'' 

The first evidence for the — > XfoK^ decay was 
reported by the Belle collaboration y, who measured 
B{B^ XcoK^) = (G-Otll ± 1-1) X 10-"* on a sam- 
ple of 31.3 X 10^ BB events. Previously CLEO had re- 
ported an upper Hmit of B{B^ XcqK^) < 4.8 x 10"'' 
at 90% C.L. (3. 

This work presents the study of the B"^ XcoK^ de- 
cay using data collected by the BABAR detector operat- 
ing at the PEP-II asymmetric energy e~^e~ collider. The 
data sample consists of 81.9 fb"^ collected at the T(4S') 
resonance containing 88.9 x 10^ BB pairs. 

The BABAR detector is fully described elsewhere . It 
consists of a tracking system for the detection of charged 
particles, a Cherenkov detector (DIRC) for particle iden- 
tification, an electromagnetic calorimeter and a detector 
for muon and identification. The tracking system in- 
cludes a 5-layer, double-sided silicon vertex tracker and a 
40- layer drift chamber filled with a mixture of helium and 
isobutane, both in a 1.5-T magnetic field supplied by a 
superconducting solenoidal magnet. The DIRC is a novel 
imaging Cherenkov detector relying on total internal re- 
flection in the radiator. The electromagnetic calorimeter 
consists of 6580 CsI(Tl) crystals. The iron flux return is 
segmented and instrumented with resistive plate cham- 
bers for muon and K'^ identification. 

Events with BB pairs are selected by requiring the 
presence of at least three charged tracks, the ratio of the 
second to the zeroth order Fox- Wolfram moment to 
be less than 0.5 and the total energy of all the charged 
and neutral particles to be greater than 4.5 GeV. We 
consider only events where at least one track identified 
as a kaon has a momentum greater than 900 MeV/c in 
the e+e~ center-of-mass frame. 

We reconstruct the XcO meson in the decay modes 



XcO ^ TT+TT" and XcO ^ K^K^ from an oppositely- 
charged pair of tracks identified as both pions or 
both kaons, respectively. Candidates for the decay 
B^ XcqK^ are formed by combining a track identified 
as a charged kaon (referred to as the "bachelor" kaon in 
the following) with a XcO candidate and performing a ge- 
ometrical vertex fit. The efficiency for the kaon selection 
used is between 70% and 90%, depending on momentum, 
while the probability for a pion to be misidentified as a 
kaon is below 5%. All the tracks are required to have 
polar angles in the region 0.35 < 9 < 2.54 rad, to have at 
least 12 hits in the drift chamber and a transverse mo- 
mentum with respect to the beam direction larger than 
100 MeV/c. In addition, tracks consistent with being 
from K'^ — > TT+TT", rj — > Tr+Tr^Tr", A pTr~ {A —> pK~^) 
decays and tracks from 7 conversions are rejected. The 
XcO candidates are required to have invariant mass in the 
range 3.32 < m^.o < 3.50GeV/c2. 

To reject the large combinatorial background coming 
from continuum qq events, a Fisher discriminant !F [8j 
is used, built from a linear combination of 11 quantities 
related to the event shape or the B kinematics. The 
coefficients are determined by maximizing the separation 
between signal and continuum background on simulated 
events. 

The selection of B candidates relies on the kine- 
matic constraints given by the T(4S') initial state. 
Two variables are defined: the beam-energy sub- 
stituted mass, rriES = \/ {s/2 + po • p)^/Eq — |pp, and 
AE — E* — where p is the momentum of the B 

candidate, (i?o,Po) is the four-momentum of the initial 
state, in the laboratory frame, and E* is the B candi- 
date energy, ^/s is the total energy, in the center-of-mass 
frame. For the B"^ XcoK^, XcO ^ K^K~' mode, when 
an ambiguity arises in cases with the same three final 
state kaons, we select as the bachelor kaon the one with 
the highest center-of-mass momentum. 

The values of the cuts for niEs and AE are de- 
termined by an optimization procedure aimed at max- 
imizing the value of S/\/ S + B. The number S of 
signal candidates and B of background events surviv- 
ing the selection are estimated on samples of simulated 
events and data from the AE "sidebands" of the toes- 
AE plane, respectively. The sidebands are defined by 
5.2 < mEs < 5.3GeV/c2, 0.1 < \ AE\ < 0.2 GeV. The rel- 
ative normalization of the signal and background sam- 
ples is determined by assuming the value measured by 
Belle for B{B^ — > XcoK^) and the world average for 
Bixco TT+TT-) and B{xcO ^ K+R-) 0. The signal 
regions in AE and toes are defined by —45 < AE < 
35MeV, TOES > 5.2750 GeV/c^ for the XcO ^ tt+tt- mode 
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and by -70 < AE < 60MeV, mes > 5.2735 GeV/c^ for 
the XcO K~^K~ mode. 

We apply a veto on fully reconstructed B D'-*^h de- 
cays, where h denotes a tt, K or p meson. We reject B 
candidates if at least one of their decay products also 
contributes to the reconstruction of a i? — > D^*^h de- 
cay, with |A£;| < 30MeV and tues > 5.27GeV/c^ To 
reduce the residual contamination from other B decays 
with charmed or charmless mesons in the final state, we 
require the invariant mass of the pair formed by the bach- 
elor kaon with the oppositely charged track from the XcO 
decay to be greater than 2 GeV/c^. 

The main source of non-combinatorial background re- 
maining after the selection described comes from non- 
resonant B decays with the same final state as the sig- 
nal, B± K^TT+TT- and B^ K^K+K'. A reliable 
evaluation of the expected contamination from these pro- 
cesses cannot be obtained based on the available mea- 
surements. These modes are expected to behave as 
"peaking background" , that is to peak in rriES and Ai5, 
while the distribution of m^^^^ is expected to be flat: this 
is used to separate their contribution from the signal by 
means of a fit to the data, as described below. 

The background from misreconstructed XcO decays to 
other modes is studied on simulated events and found 
to be negligible with respect to the other background 
sources for both the tt+tt^ and the K^K^ modes. 

The number of signal events is extracted by a simulta- 
neous unbinned maximum likelihood fit to the jties and 
m^^^o distributions for the events in the Ai? signal band. 
Three components are assumed to contribute to the se- 
lected sample: a signal component, modeled with a non- 
relativistic Breit-Wigner function convolved with a Gaus- 
sian distribution in m^^^ , and a Gaussian distribution in 
TTiEs; a combinatorial background component, modeled 
with a fiat distribution in to^co ' and an Argus threshold 
function 10] in toes! a peaking background component, 
modeled with a fiat distribution in m^^g , and a Gaussian 
distribution in tties , assuming the same resolution as for 
the signal. 

In the fit the B^ and XcO masses are fixed to their PDG 
values 01; the XcO width is fixed to the value recently 
measured by E835 [ill, r(Xco) = (9.8± l.OiO.l) MeV/c^. 
The width of the Gaussian peak in tties and the m^^Q res- 
olution are determined from Monte Carlo samples. The 
Argus shape parameter and the relative weight of the 
three components are left as free parameters in the fit. 

We verify the goodness of the fit with the three- 
component model using a Monte Carlo technique. For 
each of the two XcO decay modes, we simulate a num- 
ber of experiments by randomly generating samples of 
events distributed in tties and m^^^ according to the dis- 
tributions used in the fit. The number of events gener- 
ated for each sample is equal to the number of events 
in the corresponding real data sample; the parameters of 
the distributions are set to their fixed or fitted values. 



For each sample, the fit is repeated in the same condi- 
tions as on real data. The pulls for the number of signal 
and background events are distributed as expected. The 
probability of having a worse fit than the one on data 
is found to be about 65% and 27% for XcO — > tt+tt" and 
XcO K^K~ , respectively. 

We check the reliability of the yield extraction on a 
sample containing known amounts of combinatorial back- 
ground, peaking background and signal events. We also 
verify the stability of the fit results against variations of 
the parameters fixed in the fit by fioating them one at a 
time. 

The signal and background yields resulting from the fit 
to the data are reported in Table [J The maximum cor- 
relation we observe is about —40%, between the number 
of signal and peaking background events, for both the 
7r+7r~ and the K^K~ modes. Figure ^ shows the mES 




5.26 5.28 

lUp^ (GeV/f ") 



3.45 3.5 
^KK (GeV/c") 



FIG. 1; Extraction of the signal yield for Xco — > 7r"'"7r~ (top) 
and XcO — > K'^K~ (bottom). Left: rriEs distribution; right: 
m^^o distribution. Dots with error bars represent the data; 
lines represent projection of the fitted functions for the three 
contributions: combinatorial background (dotted), peaking 
-I- combinatorial background (dashed), signal -\- backgrounds 
(solid). 



TABLE I: Number of signal (A^sig), combinatorial background 
(A^comb) and peaking background (Apkg) events obtained by 
the fit described in the text (with statistical errors only). 



Mode 


Asig 


comb 


Apkg 


XcO tt + tt" 
XcO ^ K+K- 


33.0t^;? 


ml- 
102II1 


^^••^-6.3 

22.2l?j 
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and the m^^g distributions for events in the AE signal 
region for the two modes considered. The results of the 
fit are superimposed. 

We evaluate the systematic uncertainty to be at- 
tributed to the yield extraction by varying each fixed 
parameter, one at a time, by its error and repeating the 
fit. This results in a 2.4% (3.3%) fractional uncertainty 
for the XcO tt+tt" (xcO K^K~) mode. 

The statistical significance of the signal, defined as 
^2 log Cmax/ where Cmax/t^o is the likelihood ratio 
for the fit with respect to the null signal hypothesis, is 
8.1 (6.8) standard deviations (cr) for the XcO — * tt+tt^ 
(XcO K^K~) mode. When the systematic uncertain- 
ties are taken into account, the significances of the signals 
become 7. 7 a and 6.4cr, respectively. 

An alternative fitting method is employed to cross- 
check the results. An unbinned maximum likelihood fit 
to the m^^Q distribution only is used to extract the yield 
for the events selected in the ui-es-AE signal region. The 
signal component is modeled with a Breit-Wigner shape 
convolved with a Gaussian resolution function, while the 
background component is modeled with a linear func- 
tion. This fit yields A^sig(xcO ^ tt+tt") = i^.^t^l and 
A^sig(XcO K+K-) = 2^.7tl% Both values are com- 
patible with the results obtained with the primary fitting 
method. 

The overall selection efficiency, estimated by using sim- 
ulated data, is (27.4 ± 1.5)% for the XcO ^ tt+tt^ mode 
and (22.3 ± 1.3)% for the XcO ^ K+R- mode. The 
quoted uncertainty is mainly due to differences between 
real and simulated data in the track reconstruction ef- 
ficiency and in the efficiency of the AE cut for the 
XcO ^ TT+TT^ mode and to differences between real and 
simulated data in the track reconstruction and particle 
identification efficiencies for the XcO K^K~ mode. 

We derive the branching fractions asB — Ngig/ (eiY^i ) , 
where e denotes the overall signal efficiency and 
Ng± is the total number of mesons pro- 

duced in the data sample considered. The value 
of iV^i is determined from the measured number 
of BB pairs, iV^^ = (88.9 ± 1.0) x 10^ and using 
S(r(45) ^ B+B-) = (0.513 ±0.013) 0. We obtain: 

B{B^ Xco(7r+7r-)i<r±) = 1.32;°:^^(stat) ± 0.09(syst), 



B{B^ ^ Xco{K+K-)K^) = 1.49t^j^(stat)±0.11(syst), 

expressed in units of 10^^. The systematic error com- 
bines the uncertainties from the determination of the 
number of BB pairs, from the branching fraction for 
T(45) B^B^ , from the yield extraction and from the 
signal efficiency. 

The ratio of the branching fractions for the XcO into 
the two modes is: 

S(XcO TT+vr-) 



TABLE II: Branching fraction for ^ XcoK^ with the 
two XcO decay modes. The meaning of the quoted errors and 
of the combined measurement is explained in the text. 



Mode 


B(B± ^ 


XcoK^) (xlO-*) 


XcO 7r+7r~ 


QQ+0.60 


±0.19 ±0.52 


XcO ^ K+K- 


9 f;Q+0-63 


±0.19 ±0.43 


combined 


2.7 ±0.7 





which is compatible within the quoted errors with the 
world average 0. 

Using B{xcO TT+TT-) = (4.68 ± 0.26 ± 0.65) x 10"^ 
and B{xco ^ K+K^) = (5.68 ± 0.35 ± 0.85) x lO^^, as 
reported by the BES collaboration we measure the 
values of B{B^ XcoK^) reported in Table HH There 
the first quoted error is statistical, the second is sys- 
tematic and the third accounts for the uncertainty on 
the XcO decay branching fractions. The measurements of 
B{B^ XcoK^) obtained with the two XcO decay modes 
are compatible. Taking into account the correlated er- 
rors, we have combined them and derived the value re- 
ported in Table ini 

In summary, we have studied the process 
B^ XcqK^, reconstructing the XcO meson 
through its decay modes XcO ^ K^K^ and 
XcO ^ "'^1' ; the measured branching fraction is 
B{B^ XcoK^) = (2.7 ± 0.7) x 10"^. The result is 
significantly different from the zero value expected from 
the color-singlet current-current contribution alone. 
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0.88t^;2?(stat) ± 0.07(syst) 



Bixco ^ K+K-) 
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